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Abstract 

The main difficulty in utilizing the 0(4) symmetry of the Hydrogen atom in practical calculations 
is the dependence of the Fock stereographic projection on energy. This is due to the fact that the 
wave functions of the states with different energies are proportional to the hyper spherical harmonics 
(HSH) corresponding to different points on the hypersphere. Thus, the calculation of the matrix 
elements reduces to the problem of re-expanding HSH in terms of HSH depending on different points 
on the hypersphere. We solve this problem by applying the technique of multipole expansions for 
four-dimensional HSH. As a result, we obtain the multipole expansions whose coefficients are the 
matrix elements of the boost operator taken between hydrogen wave functions (i.e. hydrogen form- 
factors). The explicit expressions for those coefficients are derived. It is shown that the hydrogen 
matrix elements can be presented as derivatives of an elementary function. Such an operator 
representation is convenient for the derivation of recurrency relations connecting matrix elements 
between states corresponding to different values of the quantum numbers n and I. 
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I. INTRODUCTION 

About seventy ears ago, Fock in his paper Q demonstrated that the wave functions of 
the hydrogen atom in momentum representation are proportional to the four- dimensional 
hyperspherical harmonics (HSH). From the point of view of the group theory, it means that 
the hydrogen wave functions possess the 0(4) symmetry. The source of this remarkable 
property of the quantum two-body Coulomb problem lies in the presence in of an additional 
conserving vector quantity - the Laplace-Runge-Lenz operator. 

The applications of the Fock's method to various problems of quantum mechanics are 

n 

very extensive. Lieber used it to calculate the Lamb shift in the Hydrogen atom [2] . Group- 
theoretical properties of 0(4) sjTumetry of the hydrogen atom have been analyzed in jsf. 



Schwinger showed that the hydrogen Green function can be presented in a simple form 
as the series of four- dimensional hyperspherical harmonics. Fock's method has been used for 
the calculation of the retardation effects in two-photon bound-bound transitions P]. The 
four-dimensional HSH can be chosen as a Sturmian basis set in the many-center Coulomb 
problem 0, fl 3 J S H H Q, Q, Q, Q, Q, H Q- Recently^he Fock's projection 
method has been applied to the theory of anisotropic excitons, see jl9|. 

As was pointed by Lieber P|, the main difficulty in applying Fock's symmetry to the 
specific problems is the dependence of the arguments of HSH on the energy. This prohibits 
the direct application of the Wigner-Eckart theorem to the calculation of matrix elements 
between different states of the two-body Coulomb system. Often, this problem can be 
avoided by using the above mentioned set of Sturmian functions. Those functions (denoted 
sometimes as "Sturmians") have the same form as the bound-state wave functions of the 
hydrogen atom whose arguments are chosen to be equal. Thus, it is possible to apply the 
technique of the angular momentum theory to the calculation of integrals involving Sturmian 
functions which belong to the same set 

If one needs to calculate integrals between Sturmians from different sets, there appears 
the same problem as in the calculation of integrals involving hydrogen wave functions with 
different energies. Namely, the Wigner-Eckart theorem can be applied only for integrals 
containing HSH whose arguments correspond to the same point on the hypersphere. How- 
ever, hydrogen wave functions with different energies (or Sturmians from different sets) 
correspond to different points on the hypersphere. In this paper we show that this problem 
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can be solved using the multipole expansion theorems for four-dimensional HSH. Mathe- 
matically, we have derived the connection between sets of HSH corresponding to different 
stereographic projections. That is, HSH depending on some point of the hypersphere can 
be presented as a series of HSH depending on some other point which corresponds to the 
different stereographic projection of the three-dimensional space onto the four- dimensional 
hypersphere. The coefficients in those expansions are, in fact, the form-factors of the Hy- 
drogen atom. 

The form factor integrals are also known as "generalized oscillator strengths" (GOS). GOS 
play an important role in the theory of collisions of fast particles with atoms and molecules. 



Starting from the pioneering work by Bethe 21[ , there is a large number of papers addressing 
the properties of GOS of the hydrogen atom (see review 221). In the book j23| the compact 
expressions for GOS corresponding to the transitions from Is and 2s to n, / states are given. 
The hydrogenic form-factors for transitions between arbitrary bound states were presented 
in [2^ in terms of combinations of hypergeometric functions of two variables. The integrals 
between the hydrogen wave functions and the operator of the electromagnetic interaction 
have been calculated in [2^ using the set of parabolic coordinates. The expressions for the 



bound-free matrix elements for the hydrogen atom may be found, e.g. in [26, 

Q- In all 

above works the calculations have been performed by means of straightforward evaluation 
of the three-dimensional integrals in spherical or parabolic coordinates. 

Below we demonstrate that the hydrogen form-factor integrals can be explicitly calculated 
without evaluating any integrals. Our approach is based solely on the addition theorems for 
four-dimensional HSH which have been derived using the differential technique [2^. As a 
result, we can write the hydrogen form-factor integrals as a combination of multipole terms 
(they are equivalent to the radial integrals). Each multipole term can be presented in a 
compact closed form containing derivatives of a some elementary function. This differential 
representation is very convenient for analyzing various properties of the integrals such as 
their asymptotic behavior. Moreover, using the differential representations we have derived 
a number of recursion relations connecting radial integrals corresponding to different values 
of their indices. 

The paper is organized as follows. In sec |H] we consider the connection between matrix 
elements in coordinate and momentum representations as well as the general properties of 
the stereographic projection. The expression for the matrix element of the boost operator in 
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terms of HSH is derived in sec. Ill Al We demonstrate that the calculation of the form-factor 
integral reduces to the problem of the derivation of a certain multipole expansion for the 
product of HSH with some scalar function. Then, in sec. Ill Bl we use Fock's method of 
stereographic projections in order to analyze the properties of the function whose multipole 
expansion has to be derived. 

In section UTTl we calculate such multipole expansion. First, in sec. IIII A[ the multipole 
expansion is derived in the tensor form, i.e. as a tensor product of HSH. Next, in sec. IIII Bl 
we calculate the multipole coefficients explicitly. In Section UTTl the properties of the matrix 
elements are considered regarding the expressions for the multipole coefficients. In sec. IIVBI 
the dipole limit of the matrix elements is analyzed in detail. Some concluding remarks are 
given in Sec. |3 Appendix ^contains the calculation of summations arising in the course of 
computation of the multipole coefficients. In Appendix^ we derive explicit expressions for 
the multipole coefficients. Recurrency relations for the multipole coefficients are considered 
in Appendix IdI 

Below we use bold font letters like p, k to denote vectors in four-dimensional space. For 
unit vectors in that space the hat is used, e.g. x, y. Three-dimensional vectors are denoted 
as p, k, etc. 

II. STEREOGRAPHIC PROJECTIONS AND MATRIX ELEMENTS 

In this section we present the relation between hydrogen wave function in coordinate and 
momentum representations. This is done in sec. Ill Al where also the expression for the matrix 
element of the boost operator in terms of HSH is derived. Next, in sec. Ill Bl we discuss the 
general properties of the Fock's stereographic mappings corresponding to HSH entering the 
form-factor integral. 
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A. The matrix elements in the momentum space 



The hydrogen wave function in momentum representation ipnimip) is defined as the Fourier 
transform of the wave function in coordinate representation, 



ipnimip) = (^2ny/'2 j ^ ''^"'^nim{r) dr, 

i^nlm{r) = ^2^y/2 j e'^'^'lpnlmip) dp, 



(1) 



where n is the principal quantum number, /, m are the angular momentum quantum num- 
bers, / = 0, 1, . . . n — 1 and m = — /, —/ + !,.../. The explicit expression for the wave function 
in momentum space has the form Q|, 

4/55/2 

'ipnlmip) = (-^2 ^ (-p^2)2 ^"-l^"^(^)' P = Z / U. (2) 

Here, l^n-i«m(y) is the four-dimensional HSH and y is the unit 4- vector whose components 
are defined by 



Functions Ynim are normalized to the unity 

j ynlmiy)^n'l'm'{y) d^y = 5n,n'5l,l'5m,m'- (4) 

Here, it is supposed that the parameters /3 entering the definition of y do not depend on 
indices n, n' . The connection of the surface element of the four- dimensional hypersphere dVl 
with the three-dimensional volume element dp has the form 

Note that, in the case of hydrogen, the vector y depends on the rank n of HSH (see definition 
((21) of the parameter (3). The dependence of (3 on the principal quantum number n explains 
why the orthogonality relations cannot be directly applied to the calculation of matrix 
elements. 

We remark that the method of Sturmian functions has found wide use in various problems 
of the quantum chemistry and the theory of interaction of photons with atoms and molecules. 
In contrast to the set of the bound-state wave functions of the hydrogen atom, the set of 
Sturmian functions is complete. Sturmians in momentum representation are defined by 
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eq. (j2I) in which a must be a fixed number independent of n. Thus, the parameter a is the 
same for every function from the same Sturmian set. 

In many problems of the coUision theory it is necessary to calculate the form-factor 
integral, 

^(^) = / e^'^'V'n.lr) dr, (6) 

where, for the sake of shortness, we have omitted the angular momentum indices which are 
I'm! and Im for the final and initial states Uj and n,, respectively. Using above eqs. (0) the 
matrix element F{}z) can be expressed in terms of momentum space functions as 

^(^) = 77T^ \ [ dp, [ dp^rnXvi) e-^^^^-V^^-^>„,(p2) e^^^-n (7) 



Using the definition of the three-dimensional Dirac (5-function, one can evaluate integrals 
over p2 and r in closed form which yields 

F{k) = J i^nXp- k) dp. (8) 

Using explicit expressions (0) and (0) for the wave functions and introducing the indices 
n = Ui — 1 and n' = Uf — 1, the form-factor becomes 

F{k) ^ F^n',ni = / y.W(y) , / ;^% Y^U^) dn,, (9) 
p ' J [a'^ + [p — kyy 

where the parameters a and (3 have the form 

z z ^ z z 



and the components of the unit vector x are defined by 

2a{p — k) — {p — k^ 
o? + (p — ky o? + — ky 

B. The properties of the stereographic projection 



(11) 



The geometrical meaning of the unit vectors y and x in eq. and their connection to 
the three-dimensional vectors p and k is demonstrated on Fig. ^ It is seen that there is 
one-to-one correspondence between the three-dimensional space of vectors p and the unit 
four-dimensional hypersphere. Namely, every point x on the hypersphere corresponds to 
only one vector p and vice versa. This x ^ p mapping is called stereographic projection of 
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FIG. 1: The cut through the four-dimensional hypersphere. The three-dimensional p-space corre- 
sponds to the horizontal axes. The vertical zg-axis (the energy axis) is denoted as eg, the point O 
corresponds to the origin, p = 0. The point x on the hypersphere is the stereographic projection 
of the point p shifted by the vector k. 

the hyperplane (which is, in fact, the whole three-dimensional space) on the hypersphere of 
unit radius. 

As follows from eq. ©, in order to calculate the form-factor integral one has to expand 
the product 

^ Ynlmi^) (12) 



(a2 + (p _ A;)2)2 

in terms of HSH depending on the vector y. When done so, the integral F{k) calculates 
trivially using the orthogonality property of HSH. 

As a preliminary, we note that vectors y and x can be re-written in the compact form as 



x = -e + 2a^, y = -e + 2P^-, q = p-k 



(13) 



where e is the unit vector with components e = (0, 0, 0, 1) and 4-vectors p and k are defined 

by 

P = {p,P), k=ik,(3-a). (14) 
With new notations, the expression Q for the form-factor integral assumes the form 



Fn'l' 



,nl 



(15) 



/?V2 J 

Below, for simpler presentation, we will use the renormalized HSH C„(y) connected to the 
normalized HSH by Cnim{y) = (— l)"^'A/27r2/ (n -|- 1) Ynim{y)- We will also omit the projec- 
tion indices Z,m in the notation of HSH so far this will not lead to misunderstandings. 



The possibility of the efficient use of the multipole expansion technique for calculation of 
integrals stems from the surprising fact that both the products p~^C„(y) and |p — k|^^C„(x) 
satisfy the four-dimensional Laplace equation with respect to the vector p. Namely, 

Ap^C„(y) = (a,^+ ^) lc„(y) = 0, Ap^-lj^C„(x) = 0, (16) 

where the vector p has the components p = {p,Po) and the unit vectors y and x are defined 
by eqs. ()13|) in which (3 and a are constants. Note that equations ()16|) are valid for arbitrary 
values of k. 

The equations can be proved using the multipole expansions formula derived in j3| , 



= ^ (-e + 2/54) = 4 E f ) (-1)" (-T {C'n-n,(e) ® a,(p)}n. 

3 P [ P Jn P ^0 V^l/ \P J 

(17) 

Here the curly brackets mean the four-dimensional tensor products of the rank n which 
is the generalization of the conventional three-dimensional tensor product 30 1. The 
product p~"^~^C„i(p) which occurs on rhs of eq. ()17|) is the irregular solution of the Laplace 
equation. Therefore, the action of the Laplace operator Ap on eq. (fT7|) gives zero, 

^p^^".(P) = 0- (18) 

This completes the proof of the first equation in ^T^i. The second equation can be proved 
analogously noting that 

ApF(p-k) = AqF(q), 

where F is an arbitrary function and q = p — k. 

As is seen from eq. fjl6p . it is convenient to work with the product g~^C„(x) rather than 
with g~^Cn(x). In order to establish the connection between these two combinations we 
employ the auxiliary identity 

1 1 + (x ■ e) 1 



g4 2a2 q2 

which can be proved by taking the square of the second equation in ()13|) and noting that 
= x^ = 1. Applying the Clebsch-Gordan expansion technique one can prove the relation 

(X ■ e)C„,;„(x) = Cioo{x)Cnlm{^) 

= C'lOO;nim C(n-l)«"i(x) + -j= C{oo^^''/^C(„+i)Zm (x) , (19) 
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where C[qq.^J^I^ is the four-dimensional Clebsch-Gordan coefficient. At such particular values 



of parameters, the explicit expressions for those coefficients are given in ref. 
the product g"''C„(x) can be written as 



28[. As a result, 



1 



2a2 g2 



^nlm\^) "1 2(n + 1) ^{n+l)lni\^) 



2(n+l) <-(n~i)«m(x) I . (20) 



Thus, the problem reduces to the calculation of the multipole decomposition of the function 
q'~^Cn(x). Once such the decomposition is known, the expression for the function g~'^C„(x) 
can be derived using eq. (pUj) . 



III. THE MULTIPOLE EXPANSION OF TRANSLATED HARMONICS 

In this section we derive the expansion of the function q^^Cnix) in terms of HSH depend- 
ing on the unit vector y. In sec. IIII Al the tensor form of the multipole expansion is derived 
by consequent application of the multipole series given in [2^ . In sec. IIIIBI the coefficients 
of that expansion are calculated in an explicit form. 



A. The multipole expansion in terms of tensor products 

As was noted above, it is sufficient to calculate the expansion of the product q~'^Cn{^)- 
The calculations are performed in three steps. First, we expand the product g~^C„(x) in 
terms of HSH depending on the hyper-angles of the vector q. Second, the resulting HSH 
is expanded in terms of HSH whose arguments are the hyper-angles of the 4D-momentum 
vector p. The final, third step, is to expand HSH depending on p in terms of the HSH 
depending on the unit vector y. Clearly, at each step the series of irreducible tensor products 
of HSH will occur. However, as we will see below, only one of the three series is infinite. 
Further, the Clebsch-Gordan coefficients inherent to the irreducible tensor products arising 
during the computations have simple structure and can be written in a closed form. 

We begin the calculations by using the expansion similar to eq. ()17j] . 

^;a(x) = ^y^f 'M(-l)™f-^" 



I y ni=0 ^ 1^ ^ ^ ^ 



Now, we have to re-expand the function q ^C'n^(q) in terms of HSH depending on the 
vector p. The use of the muhipole expansion formula (eq. (65) of ref. 2^) yields, 

^^-(^) = E ^^("' ^2 ^ 0a....(p)}„. (22) 

y n2=0^ ^ 2 / 

Third, we have to go back from vectors p to vectors y- In order to do that we note that the 
terms containing vector p can be re-written as 

1 2("i+"2)/2 ^ w r n 1 

^^-+-(P) = >^{P}-+- = 2^^^^""^^/^ ^ • (23) 

f t' (.f J ni+n2 

The connection (fT^ between vectors p and y leads to the chain of equations 
2{ni+n2)/2 f P_l ^ I ie + vl ^ 

^ ni+n2 / + n \ 

= (2/?)»i+"2 E (^'''^^'''j{Cni+n2-n3(e)®C„3(y )}„,+„,. (24) 

Finally, we have to substitute the above equations (j22I), and into eq. (PT|) . As a 
result, we arrive at the identity 

ic„„„(x) = (ziT ^ ( " ) ,_„)... + + 1) "f " + 

^ ^ „i=0 ^ n2=0 ^ ^ / „3=o ^ 

X {C„_„^(e) O {C„2(k) (g) {C„,+„2_„3(e) ® C„3(y)}„,+„2 (25) 
where the short-hand notations have been introduced 



u=-, v = — = . (26) 

The above equation gives the connection between the hyperspherical harmonics C„(x) 
and C„3(y)- As is seen, among three summations only one (over the index 712) is infinite. 
Expansion is the farthest point one can reach using the angular momentum technique. 

If we substitute the expansion (j2SI) into the integral (|TH|l . then, due to the orthogonality 
of HSH, only the terms with = n will remain. Thus, the matrix element will be defined 
by the two summations one of which is infinite. Below we demonstrate that the infinite 
summation can be evaluated in closed form as the derivatives of a some elementary function. 
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B. Explicit expressions for the coefficients of the multipole expansion 



The tensor product in the expansion ()25|) has a remarkable property. Namely, the internal, 
as well as external, tensor products are "minimal", i.e. the rank of the product is equal 
to the sum of ranks of HSH entering the product. The intermediate product having the 
rank ni, contains tensor multipliers with the ranks n2 and ni + n2. Such tensor product 
can be considered as "maximal" since its rank is equal to the difference of the ranks of 
the constituent tensors. It turns out that the four- dimensional Clebsch-Gordan coefficients 
corresponding to the minimal or maximal tensor products can be evaluated in the closed 
form, see 28[. As a consequence, the tensor construction in eq. (I25j) can be significantly 
simplified. 

The calculation of the three-fold tensor product in eq. (j^^ simplifies also by the fact that 
the hyperspherical harmonics of the unit basis vector e reduces to a constant number 



CnZm(e) = Vn+ I5ifi5mfi- (27) 

Using this equation and the definition of the tensor product of HSH j3| , we can write the 
tensor product from the right-hand side of eq. as 



'^nlm ^ = {C„_„^(e) ® {C„2(k) ® {C„j4.„2-ra3(e) ® C*ra3(y)}ni+n2}ni}reZm 

ilm ^{ni+n2)hms 
2(2012; (ni+»i2)'3"i3 (rai+n2— ".3)00; ri3'3'Ti3 



E\ ^ f^nlm ^riilm ^(ni+n2)/3ni3 

/ ^ '-'(n-ni)OO; niim ^n2l2m2; (ni+n2)i3m3 C 



l2,h "12,^3 



X A/(n -ni + 1) (ni + 7^2 -n^ + l)Cn2l2m2{^)Cnzhm3{y), (28) 



where Cl^'iJln'-n"i"m" denotes the Clebsch-Gordan coefficients of the 0(4) group. 

We have to calculate the summation over the index n2 in (j2Hl). First, we note the explicit 
expression for HSH C„2i2m2(k), 



Cn2l2m2i^) = HY^f^^^f {260)0,2^2(0 A), (29) 

where Ci^rn2{0i<P) are conventional spherical harmonics depending on the spherical angles 
6*, of the vector k. The hyper- angle Oq is defined by 

a ^ ~ " . a |fc| 

coso'o = — - — = — , , smfc'o = — , (oUj 
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The functions Xi'^^'^ in eq. ()29|) are the generahzed characters of the 0(3) rotation group 



which are connected to the Gegenbauer polynomials 
HSH Cn2i2m2 has the form, 



. Thus, the explicit expression for 




where C^^i^ is the Gegenbauer polynomial. Substituting this equation and eq. into 
eq. (P^jl we arrive at the multipole expansion 

I (—1)" °° ^ ^ 

"3=0/3=0 Z2 m2,m3 

where Cj^^i^^^ is the three-dimensional Clebsch-Gordan coefficient, and the index I2 runs 
from |/ — /3I to / + /a so that the sum / + /2 + ^3 is an even number. This is due to the fact 
that the coefficients C^q^^q vanish at odd values of the sum I + I2 + h- The coefficients A in 
the above equation are defined by 



where the summations are performed over all non-negative values of indices rii and n2 so 
that all binomial coefficients and factorials remain finite. As is seen, the summation over 
rii is finite while the summation over n2 is infinite. These summations are calculated in 
Appendix ^ The sum over I2 in (jH^ has simple physical meaning. It is the multipole 
expansion with respect to the "momentum transfer" vector k. 

As was pointed above, the multipole expansion for the product g~'^C„(x) can be derived 
from the expansion for the function g~^C„(x) by applying eq. ((201) • Omitting details of some 
routine transformations, we present the final result, 

I ( — 1)" °° "3 I (i ) 10 



^ ^ 713=0 ^3=0 I2 



where the coefficients B^]^. 



y^ Cl^2l3m3^l2m2{(^A)Cnsl3ma{y), (34) 



X 

"12, "^3 




,(/2) ^_(^\ ^2' l^t2 ^ / (n + /+l)!(n3 + /3 + l)! (;^) 



12 



and the functions P^f^s h defined by eq. ()A9|) of Appendix |X1 

At zero "momentum transfer" (i.e. at k = 0), the properties of the three-dimensional 
Clebsch-Gordan coefficients in the expansion (jH^ lead to the vanishing of all terms with 
I2 > and I ^ h- Thus, the expansion (jH^ reduces to 

1 °° 

-iCnU^) = ^ E(-l)"^"^<;n3,z^n3.„.(y), (36) 

713=/ 

where the parameter B^^\.^,^ ^{k = 0) is calculated in Appendix iBl 

For the sake of completeness we note that the multipole expansion ()25|) can also be written 
as the series containing the four-dimensional sohd harmonics, Cn2i2m2^ = ^"^C*n2«2m2(^)' 

— Cnlm{^)= 2 C 10 120^ l2rn2 hms^ nahma (y) Cn2l2m2 0^) Dn2ni, 113(0:, P), (37) 

where the coefficients D„2n3,/,i3(«, /?) are defined by 

n m (ni + n2 + /3 + l)!(^i+^2-/3)! /mm 



ni 



D=r-lVV2m-- / (n + / + l)!(n-/)!(2/, + l) 

^ ' ^ ^' (r^3 + /3 + 1)! (r^3 - /a)! (^2 + ^2 + 1)! (n^ - h)\ 



(38) 



The summation over rii in eq. (jHHj) can be calculated in closed form which yields the hy- 
pergeometric function 3F2, see ref. j^^. Thus, the expression for the coefficients Dn2n-i,ih 
becomes 

[n2 + l + h + l)\{n2 + l-h)\ 



D^2n3,iUo:, f3) D _ + 1), 

/^\' / /-n, n2 + /-/3 + l, n2 + / + /3 + 2; -f \ 

I • (.•^yj 

^2 - ^3 + / + 1, 2/ + 2 

The analysis of the series ()37|) shows that the region of convergence with respect to the index 
n2 is defined by the two conditions, \a\ < and \{a — [3Y + k'^\ < 2\[3\. The series in (P7|) 
diverges when one of parameters a, (3 is purely imaginary and another one is real. 

The convergence of the above expansions (j34|) and (j37jl with respect to the index 77,3 of 
HSH is practically not important because these expansions in will always be truncated 
after the integration over the hyper-angles of the unit vector y. 
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IV. PROPERTIES OF THE MATRIX ELEMENTS 



In this section we derive the exphcit expressions for the matrix element in eq. (jlSp in 
terms of the coefficients B^f.^^ and P^fna I3 introduced in the previous section. In sec. II V A) 
the general formulas are given and few examples are considered. In sec. II V BI we analyze the 
long-wave (i.e. dipole) limit of the matrix elements. 



A. Matrix elements in general case 

In order to derive an expression for the form-factor matrix element, we replace in the 
decomposition HSH Cnim with the normalized HSH Ynim and insert the resulting ex- 
pression into (|T3|). This leads to the equation 



Fn'l',nl — 

n3,l3,h 



X E ClTm,i3m3Ci,m,ie,<l>) j , ^, {y)Y^,i,^,{y) d^ly . (40) 

m2,ni3 

Since HSH are orthonormal, the integral over the hypersphere is equal to unity for n^l^m^ = 
n'l'm' and is zero otherwise. As a consequence, we obtain the expression 

Fn'l'M = \ ni 1 ^ ^'''^ro°2oQ2'^2/'m''^'2m2(6',0)5fjJ^,;,. (41) 

V ' h=\l-l'\ 

If the transitions between the states of the hydrogen atom are considered then the factor 
under the square root reduces to unity as follows from the definition (|10|) of parameters a 
and p. 

The above formula for the form-factors Fniu,ni can be re-written in terms of radial inte- 
grals. Using explicit expressions for the radial hydrogen wave functions and decomposing 
the plane-wave in eq. (jH)) over the three-dimensional spherical functions, eq. ()41|) can be 
brought to the form 

^ ^^+''+2e-(«+/3)r-^-^^ j^l^l^j $(^/ _ ^/^ 21' + 2; 2(3r) $(/ - n, 2/ + 2; 2ar) dr 

(42) 



,,+,^+,+1 y(2/ + l)!(2r + l)! / \k\\'" ^(,,) 



(2a)'+i(2/3)''+2 \ p 
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where the parameter w is defined by eq. f)A10|) of Appendix $(/ — n,2l + 2; 2ar) is the 
confluent hypergeometric function and ji^ is the spherical Bessel function. The definition 
(jA9|) of functions P^fn' u imphes the following limitations on the indices in eq. (fl^ , 

l<n, I' <n', l-k + l' + 1>0. (43) 

(The requirement for the combination I + I2 + I' to be an even number is not necessary 
in eq. (02)) •) Thus, we have shown that the integrals of the kind fl42|l can be written as a 
three-fold derivative of an elementary function (cf. eqs. (jl^ and ()A9jl ). It turns out that 
for small values of indices n and / (or n' and /') the functions P^^Jn' v ^^"^ written in a 
compact form as a combination of Gegenbauer polynomials. 

For example, the matrix elements for the transitions from the ground state n = / = can 
be written as 

I 



where Rn^i are normalized radial wave functions and the function Pq g. ^ ^ is the combination 
of two Gegenbauer polynomials (see eq. ()C6|) of Appendix El) • We remind that the index n 
above is equal to the principal quantum number minus unity. For hydrogen matrix elements 
one has a = 1 and (3 = l/(n + l) while for Sturmian matrix elements the parameters a and (3 
can be arbitrary numbers. As an example, on fig. |21the parameters |-Fn,«P corresponding to 
transition Is ^ 5/ are plotted as functions of In (The reason for such choice of abscissa 



variable is explained in [2^ |33|.) As is seen, the hydrogen matrix elements have strong 
maxima at some values of \k\. The Sturmian matrix elements behave in less regular manner, 
e.g. they can have several peaks at different values of the momentum \k\. 



B. The dipole limit 

In the dipole limit we have that — > and, hence, it is necessary to retain only the terms 
of zeroth- and first-order with respect to Such terms correspond to the two values of 
the multipole index, ^2 = and I2 = l- The derivation of the parameters P^^ini 11 for such 
values of I2 is considered in Appendix O 

According to eq. ()C12|) of Appendix [01 the parameters with I2 = can be written in 
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FIG. 2: The parameters |-F4,ip. (a) the matrix elements for the hydrogen atom: a = 1, /? = 1/5. 
(b) the matrix elements for Sturmians: a = /3 = 1. 



closed form in terms of Gauss hypergeometric function 

(0) ^ 4(-l)"' / -An \ (l-u\ n'+l-{n + l)u 

n,i;n',«W {21 + l)\ \{l - ) \l + U ) {I + uf 



-Au 



X 2F1 [l-n,l- n'- 21 + 2; (^3^J • (45) 
The Gauss function in this equation is a polynomial whose order is min(n — l,n' — I). It is 

(0) 



easy to see that at m = {n' + l)/(n + 1) (as well as at u = 1) the parameter -P„/„/ ; is 



and the form-factor matrix element becomes 



Fn'l' 



n'l'nl 



k=0 



(n-iy. 

(n + l + l) 



(46) 



(47) 



This reflects the orthogonality property of the hydrogen wave functions (or, at m = 1, 
Sturmian functions from the same set). 

The parameter P^ln' i-i describes the dipole transitions and expresses as (see eq. ()(^13|1 
below) 



2f-l)" 



A \ l+l / T \ n+n' 



(n' + / + l)(n' + / + 2)/o 



3(2/ + !)! V(1-m)V \1 + u, 

+ 2{n' + I + l){n' - I + + {n'-l + l){n' - 1)^ k (48) 



where the coefficients fk are defined by 

2~k 

I I — II. \ 

fk 



l + U 



{n + l)u-n' -2 + k ( , -Au \ , , 
^ ^ 2^1 ( / - / - 1 - + fc; 2/ + 2; — ) . (49) 



(1 + m)2 



{\-uf 
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In the case of matrix elements between the hydrogen states, we have that u = (n' + l)/(n + l) 
and the coefficient /i vanishes. Thus, the matrix element in this case is a combination of 
two Gauss hypergeometric functions in agreement with the well-known Gordon formula js^ . 
The expression for the coefficient -P^';.^„/ can be derived from eq. (j48|) using the symmetry 
property ()A11|) which in the dipole case has the form 

^il',.W = «"^"'^^'^^^S>/l/«)- (50) 

V. CONCLUSION 

In the presented paper we have applied the method of Fock's stereographic projections 
to the calculation of matrix elements involving the hydrogen-type wave functions. We have 
shown (sec. |n} that the problem of the calculation of matrix elements is equivalent to the 
derivation of a multipole expansion of the four-dimensional hyperspherical harmonics with 
their arguments defined by means of the shifted stereographic projection. 

There is a conceptual similarity between our approach and the work 35]. In that paper 
the two-photon transitions between the states of the hydrogen atom have been considered 
using the representation of the Coulomb-Green function with two free parameters. Our 
multipole expansions of the hydrogen wave function in momentum space (see eqs. (I32|l 
and (jnH)) also contain free parameters f3 and k. By "tuning" the energy parameter /?, we 
can move the argument of HSH to any desirable point on the hypersphere. Usually, the 
choice of the parameter f3 is fixed by demanding the orthogonality of HSH entering the 



integral. The same approach was used in 



35l |. The proper choice of free parameters can 



result in the significant simplification of expressions for the compound matrix elements or 
in the imporoving the accuracy of numerical computations. The method developed in the 



3, 



36l | in order to account the 



presented paper can also be used to generalize the results of 
non-dipole effects. 

The expression for the form-factor matrix element in terms of the multipole functions 
Pnfn' I' is given in sec. IIVI We showed that the functions P^t^n' v icici^ integrals 

between the radial hydrogen wave functions and the spherical Bessel function arising from 
the multipole expansion of the boost operator exp(ifc ■ r). 

It turns out, that the closed expression for P^f^' v terms of special function of a 
single variable does not exist in general case of arbitrary values of the indices n,l,n',l',l2 
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and arguments a,P and k. Only in some particular situations, e.g. in the long-wave (i.e. 



ih) 



dipole) limit as well as for the transitions from the ground Ig-state, the functions P^l' 



can be written as a combination of Gauss hypergeometric functions Surprisingly, there 
exists a remarkably simple differential representation for coefficients P^fn' v given by 
eq. ()A9|) of Appendix ^ This differential formula can serve as a source for many explicit 
expressions for functions P^fn' v '^^^ differential representation can also be effectively used 
for the derivations of the recurrence relations for P-functions. A number of such relations 
is derived in Appendix O We note also that recursions involving radial hydrogen matrix 
elements of the powers have been considered in j^. 

The results of the presented paper can also be generalized on the case of bound-free 
transitions. In general, this can be achieved by applying the replacements n ±i/\/E, 
where E > is the energy of the continuum state. At this stage the differential operators 
in eq. ()A9|) become undefined. However, this problem can be avoided by using the contour 
representations for the derivatives. 

The matrix elements of the operator of electromagnetic interaction e exp{ik ■ r) can also 
be derived using the formalism developed in this paper. In this case one should replace 
with (e-x)g~^ in eq. (fT3j) where e = (0, e) is the four-dimensional photon polarization vector. 
Then, one can apply the Wigner-Eckart theorem to the ensuing integral. The corresponding 
results can be written again in terms of the functions P^^J^' v- 

One of the possible applications of the technique developed above could be the calcu- 
lation of the two-center Coulomb integrals. Namely, by applying the Fourier transform to 
the function l/|ri — one arrives at the form-factor integrals for which the differential 
representations can be used. The work on this problem is in progress. 
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APPENDIX A: THE CALCULATION OF SUMMATIONS IN EQ. 



In this Appendix we calculate the two summations in ()33p. We begin by calculating the 
sum over n2 which we denote as Sn^ , 

{ni+n2 + h + l)\{ni+n2-h)\ 



n2=l2 



COS 610) 



(ni + 712 - ''^3)! (^2 + h + 1)! 
Replacing the summation variable hj k = n2 — h we obtain 



(Al) 



Sn, = V'' 



t^,' {k + n,-n, + h)\{k + 2h + l)\ ■ 



{k + ni-n^ + h)\{k + 2h + l)\ 
Factorials in this identity can be removed by introducing the differential operators 
{k + ni + l2 + h + 1)! {k + ni + h- h)\ 



(A2) 



{k + ni-n^ + l2)\{k + 2l2 + l)\ 



-IT 



-{k+2l2+2) 



(A3) 



t=l,T=0 

where dt = d/dt and dr = d/dr. Inserting equation ()A3|) into ()A2|) one arrives at the series 
which can be evaluated using the generating function for Gegenbauer polynomials 

' ' 2t;cos^o 



k=0 



t + T it + T)\ 

[v^ -2vcos9Q{t + T) + {t + TY]-^^-\ (A4) 



Thus, the parameter can be written as 



r 



,21-/2-1 



Now we have to evaluate the sum over rii in eq. (jHSJ, 



n — / 



n—l 



— riij {rii + / + 1) 
11,1 — (, 

Using eq. ()A5|) we write 5"^ in its full form 



ni=0 



ni+/ 



(ni + 2/ + l)!'^"^ 



t=l,T=0 

(A5) 



(A6) 



n—l 



n — I 



u 



n\+l 



rii J {m + 21 + 1)1 



- 2v{t + r) cos 9o + {t + t) ] 



21-/2-1 



(A7) 



t=l,T=0 
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After some rearrangements using the chain differentiation rule, this identity can be written 
in closed form, 



Sri 



I2 



in + l + 1] 



X 



[v'^ - 2v{t + r) cos eQ + (t + rf] 



21 -'2-1 



■ (A8) 



t=l,r=0 



As is seen, the two-fold summation in ()33|) reduces to the derivatives of an elementary 
function. 

According to eq. (j^n|) . in order to calculate the coefficients in the expansion (jH^ of the 
function g~'^C„(x), one has to evaluate the combination of three functions Sn- Omitting 
the details of some routine transformations, we present the final result for the parameters 



n,l; n3,l3 



in dngi. 



{n + l + 1)! (77,3 + k + 1)! 

X [w' + (1 + «)(t + r) + (t + r)2] 



where the short-hand notations have been introduced 



w 



--^/l — 2v cos 60 + f ^ 



2P 



We note also the following symmetry property of functions P^t-n i ('^' 



(A9) 



t,T=0 



(AlO) 



(All) 



This equation can be proved by making the replacements t ^ ur and t ^ ut'ui the definition 



APPENDIX B: THE EXPLICIT EXPRESSIONS FOR THE FUNCTIONS , 

n,l; 713,13 

In order to calculate the functions P^^^n^ in eq. ()A9j) in an explicit form, we use the 



generating function for the Gegenbauer polynomials 

00 



+1 



9=0 



U+l 

2w 



fBll 
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Acting by the operator di '^+'3+1 on this equation one observes that the terms with q < 
I — I2 + h + 1 vanish. Thus, the functions , can be written as a combination of 

Gegenbauer polynomials 

where q'max = n + 713 — / — and parameters are polynomials in m, 

C fn) = V (^"^ ~ ~ (/-/2 + /3 + g + l)! .33. 

where the summation over m is taken over all non-negative integer values at which factorials 
remain finite. The parameters Cq{u) can also be written in the following form 

^^^"^-^ ^' (n + / + l)!(n3 + /3 + l)!g! ^* ^^^)^"3-/3 W 



t=0 



where L^n'^^^ is the Laguerre polynomial 38 1. This equation can be proved by straightfor- 



ward calculations. Thus, the parameter Cq{u) is proportional to the coefficient at the g-th 
power of the product of two Laguerre polynomials. 

One can obtain an explicit expression for Cq{u) in terms of generalized hypergeometric 
function 3F2. There are two cases which must be considered separately. For q < — h the 
index m in ()B3|) runs from to min(g, n — I). This leads to the formula 



q J{2h + q+l)\(2l + l)l. ' y 21 + 2, n,-k-q+l 

For q > — one has that g — ns + ^3 < m < n — /. In this case the above formula is 
invalid and the correct expression has the form. 



q-n3 + hj i2l-n3 + h + q + 1)! (ng + h + I) 

I + l3-n-n3 + q, h-ns, -/g - rig - 1; -u 
X 3i^2 • (B5) 

2/ + 2, h-n^ + q + l 



APPENDIX C: FUNCTIONS P^''^ , IN PARTICULAR CASES 

In this appendix we analyze the situations when the functions P^f^^^ can be written in 
closed form. Our consideration is based on the general differential relation 

at 



a — f5t 



(CI) 
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Here, a,/3,7 are arbitrary numbers and f{t) is an arbitrary function. This formula can be 
proved by expanding the function / into Taylor series and comparing coefficients at the 
expansion terms. 

Let us consider the functions P corresponding to the transition from the ground Is-state. 
In this case we have n = I = and, hence, l2 = l?j = I. According to the definition eq. ()A9|) . 
we can write 



In order to calculate the action of the differential operators we note the formula 

(1 + t)"+'+i dt = a(i + - + / + i)(i + (C3) 

Now the parameter -Pg g-n i '^^^ re-written as 

H%n,i = (^T^^Vl)! ^^"''^'^^ +t)-{n + l + l) dn (1 + + + l)t + tV-\_^- 

(C4) 

By applying the differential transformation ()C1|) to this equation we arrive at the equation, 



p!^Li = (^^TtVt)! -(^ + ^ + 1) dr') (1 - t)[w' - 2at + MV^^, (C5) 

where a = — (m + l)/2 and the parameter v is defined in eq. fl26p . The calculation of 
derivatives can be performed by expanding the term in square brackets similarly to eq. (jBljl . 
After some simple algebraic manipulations we obtain 



+ („ + ; + l)!^:!^cW_,W|. (C6) 



where the argument of the Gegenbauer polynomial is 

_ w^-iu + l)/2 _ a^-P^ + P 
vw Ap'^vw 

Similarly, one can derive rather compact expressions for the P-coefficients for the transitions 
from 2s and 2p states. However, for the sake of shortness, we do not present here the 
corresponding equations. 
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The dipole limit 

Now, let us consider the functions P^i-ns h limit of small \k\, i.e. the dipole limit. In 

this case only the zeroth- and first-order terms (with respect to |A;|) should be taken into 
account. Thus, in the dipole limit one has that w = {u + l)/2 and it is necessary to retain 
only the terms with ^2 = and ^2 = 1 in eq. ()A9|) . 

Let us calculate first the function P^lns v Since A; = 0, the term in square brackets in 
eq. ()A9|) reduces to [w + t + r]~^. This simplifies the calculations, so that eq. ()A9|) becomes 



(C8) 



t,T=0 



It is convenient to present the last multiplicand in this equation as (w + t + r) ^ 

-{l/{2l + 2))du,{'w + t + T)~^'-\ 



t,T=0 

(C9) 



Using the auxiliary identity ()C1|1 we can calculate the derivative over t in closed form, 

d^'-\w + r + t)-'^-\l + tr+'+' _ = ^""^ I ^ ^ (w + r)-"3-^-2(^ _ 1 + rr^'i, (ClO) 

t — ( I X 1 . 



By applying eq. ()C1|) to the calculation of the derivative over r one arrives at the Gauss 
hypergeometric function, 

Pi'Ui = - ir-";^'-%Fi(/ - n,/ - n3;2/ + 2;;^), (Cll) 

where z = u/[{w — u){w — 1)]. 

The calculation of the derivative over w simplifies by the fact that dw z{w) = at the 
point w = (m + l)/2. Thus, the parameter Pn}-nz,i assumes the form 

P^L,i = J^[:^_ ^ni,+3 [ns-n + 2{n+l)il-w)]z^-%F,{l-n,l-ns;2l + 2;z). 

(C12) 

In the case of /2 = 1 there are two possibilities. Since the combination I + I2 + I3 must be an 
even number, we have that = /±1. Below we calculate the coefficient P^i^^ i-ii"^^ '^)- (The 
coefficient Pni-i-nsi '^^^ ^e derived using the symmetry relation (jXTTJ.) The calculation 
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procedure is similar to that considered above. As a result, we have, 

/IN 7,"3+2 

-^n,i;n3,«-l - g (^2/ + 1)! ^ ^ * 

X (Z - n, / - 1 - rig; 2/ + 2; z) , (C13) 

t=i 

where z = ut /[{w — u){w — t)]. Omitting details of somewhat lengthy calculations we present 
only the final result which is given by eq. PHj) of sec. IIVI 



APPENDIX D: RECURRENCE RELATIONS FOR THE FUNCTIONS P^''^ , 



The differential representation ()A9|1 can be used for the derivation of the recurrence 
relations for the functions P^ins i^- Such relations can be derived using the properties of 
the differential operators. For example, we note the following operator identity 

= (ns + h + l) dt'~''~\l + wtr+'' + dt'''''\l + wtr+''+' dt . (Dl) 

(It can be proved by straightforward calculations.) A similar identity can be written for the 
operators dr- These two identities result into a pair of recurrences for the functions Pn^Jns I3 
with fixed multipole index I2, 

{^3 + ^3 + 1) Pn,l;n3,l3 ~ {^3 ~ h) Pn,l;n3-l,h — Pn,l;n3,h-l — Pn,l;ni-l,h-l^ 
{n -|- / -|- 1) Pn,l;n3,l3 ('^ Pn— 1,1, 113,13 Pn,l—l;n3,l3 Pn~l,l—l;n3,l3y 

where we have omitted the superscript I2 for the sake of shortness. The recursion connecting 
functions with different indices I2 can be derived using the differential identity, 

gi-i,+i3+2 ^ 2(/2 + l)[{i-t- t) r^^^^ _ (/ _ /2 + + 1) d[-''^'' Ri,+i] , (D3) 

where i?/2 = [1 — 2^(t+r) + (t+r)^]~'2-i and ^ = — (1+m)/(2w). After some transformations, 
the above identity leads to the recursion connecting functions P with indices I2, ^2 + 1, 



{n + l + 2) Pil{^^ + (n - /) pS,„3 - 2(n + l)Pi'^^^ = (k + 1) ((n + / + 2)(1 - u) Pj^ 



ih+l) 

+I,n3 



+ {n-l){l + u) Ptt'i + 2[u{ns + k + 2) - n - 1] P^^f^ - 2{n, + h + 2)u P^t^O . 



(D4) 
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where -P^%3 = -P„ ^ „3 , etc. The recursion connecting P-functions with fixed values of an- 
gular momentum indices l,l2,h can be derived using the fact that dr Rh — dt Rh- Omitting 
detail of the routine transformations, we present only the final result, 

(ns + ^3 + 2) 

+ 2[n+l-«(n3 + l)]Pn,n3 = 0, (D5) 

where P — P^'^^ 

Finally, we note that there are two other recurrence relations which connect functions 
P with different values of indices n, 713. Both of those recursions contain eight terms and 
they form a complete set with respect to the indices n and ^3. It means that having 10 
"neighbor" parameters one can recursively calculate the parameters Pn.ns with arbitrary n 
and 77-3. However, since these recursions have quite cumbersome form, we do not present 
them here. 
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